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ABSTRACT 



The ACS Survey of Galactic Globular Clusters, an HST Treasury Project, 
will deliver high quality, homogeneous photometry of 65 globular clusters. This 
paper introduces a new collection of stellar evolution tracks and isochrones suit- 
able for analyzing the ACS Survey data. Stellar evolution models were computed 
at [Fe/H]= -2.5, -2.0, -1.5, -1.0, -0.5, and 0; [a/Fe] = -0.2, 0, 0.2, 0.4, 0.6, and 0.8; 
and three initial He abundances for masses from 0.1 to 1.8 M Q and ages from 
2 to 15 Gyr. Each isochrone spans a wide range in luminosity from My ~14 



1 Current Address: Astronomical Observatory, Volgina 7, 11160 Belgrade, Serbia 

2 Computational Research Division, Lawrence Berkeley National Laboratory, MS 50F-1650, 1 Cyclotron 
Rd, Berkeley, CA 94720 



- 2 - 



up to the tip of the red giant branch. These are complemented by a set of He- 
burning tracks that extend from the zero age horizontal branch to the onset of 
thermal pulsations on the asymptotic giant branch. In addition, a set of com- 
puter programs are provided that make it possible to interpolate the isochrones 
in [Fe/H], generate luminosity functions from the isochrones, and create syn- 
thetic horizontal branch models. The tracks and isochrones have been converted 
to the observational plane with two different color-T cff transformations, one syn- 
thetic and one semi-empirical, in ground-based B, V, and I, and F606W and 
F814W for both ACS-WFC and WFPC2 systems. All models and programs pre- 
sented in this paper are available from http://stellar.dartmouth.edu/~models/ 
and the Multimission Archive at the Space Telescope Science Institute (MAST; 
http:/ /archive. stsci.edu). 

Subject headings: globular clusters: general — stars: evolution 



1. Introduction 

The ACS Survey of Galactic Globular Clusters (Sarajedini et al. 2007), (Anderson et 
al. 2007, in preparation) will set photometric standards for years to come. The survey has 
obtained high quality photometry of 65 globular clusters from the horizontal branch down 
the main sequence to ~0.2 M at the least. The data from this survey will allow for revised 
absolute and differential age analyses, studies of the stellar luminosity and mass functions 
and dynamical evolution of a large number of clusters in a homogeneous data set among 
many other things. The data will provide a critical test for some aspects of current stellar 
evolution calculations, in particular the shape and luminosity function of the lower main 
sequence. The survey will eventually combine ACS-WFC and WFPC2 photometry in the 
F606W and F814W bands and ground-based B, V, and I photometry for most of the clusters. 
In order to analyze these data, it is preferable to have theoretical models transformed to the 
appropriate photometric system that can be applied directly to the observations. Sirianni et 
al. (2005) provide a detailed analysis of all aspects of the ACS, including filter transmission 
curves that account for the performance of the system from optics to electronics, thereby 
allowing accurate synthetic photometry to be produced from synthetic fluxes (Bedin et al. 
2005, e.g.). 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555, under programs 
GO-10775 (PI: Sarajedini). 
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There are a number of groups currently producing isochrones similar to those introduced 
in this paper, chief among these are Berbusch & Vandenberg (2001); Vandenberg et al. (2006, 
hereafter Victoria- Regina) , Girardi et al. (2002); Bonatto et al. (2004); Girardi et al. (2004, 
hereafter Padova), Pietrinferni et al. (2004, 2006, hereafter BaSTI) and Cordier et al. (2007), 
and Yi et al. (2003, 2004); Demarque et al. (2004, hereafter Y 2 (Yonsei-Yale)). Among these, 
all provide isochrones for ground based B, V, and I bands but only the Padova and BaSTI 
libraries include isochrones on the HST photometric systems. Of these, only the Padova 
library extends down past the lower mass limit of the ACS survey data; the other terminate 
at 0.4 M (Y 2 ) or 0.5 M (Victoria-Regina and BaSTI). 

In order to analyze the full complement of photometry a new set of stellar evolution 
models have been produced. This paper describes in detail the stellar evolution code used 
to create the models in §2; the color-T cff transformations used to convert the models to 
the observational plane in §3; properties of the model database, including descriptions of 
the tracks, isochrones, luminosity functions, and synthetic horizontal branch models, in §4; 
comparisons of the isochrones with other sets isochrones and data from the ACS survey in 
§5; and finally the results are summarized in §6. Following the main text is an appendix that 
describes in detail the stellar evolution track, isochrone, luminosity function, and synthetic 
horizontal branch model files and outlines the use of computer codes distributed with the 
models. 



2. The Dartmouth Stellar Evolution Program 

The stellar evolution code described here is the Dartmouth Stellar Evolution Program 
(DSEP). Many details of DSEP have been previously described by Chaboyer et al. (2001); 
Bjork & Chaboyer (2006). This section is intended to explain the physics and numerical 
methods used in the code that have not been covered by Chaboyer et al. (2001); Bjork & 
Chaboyer (2006). 

Stellar evolution models have been created for [Fe/H]=-2.5, -2.0, -1.5, -1.0, -0.5, and 
with [a/Fe]=-0.2, 0, 0.2, 0.4, 0.6, and 0.8. In the models, a-enhancement refers to en- 
hancements in the following a-capture elements: O, Ne, Mg, Si, S, Ca, and Ti by the 
same amount as specified in [a/Fe]. The models assume the initial He mass fraction follows 
Y=0. 245+1. 54 Z where Y=0.245 is the primordial He abundance from Big Bang nucleosyn- 
thesis (Spergel et al. 2003) and the required initial He abundance of the calibrated solar 
model is Y=0.274. In addition, sets of models were produced with Yj n ^=0.33 and 0.4 for the 
same range of [Fe/H] with [a/Fe]=0 and 0.4. Complete details of the initial compositions 
for all models are listed in Table 2. 
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For future reference, the solar abundances were adopted from Grevesse & Sauval (1998) 
and global solar parameters luminosity, radius, surface composition, and age from Bahcall et 
al. (2005). Table 1 defines the requirements for a calibrated solar model and gives the values 
obtained from the DSEP solar-calibrated model. The calibrated solar model has a mixing 
length parameter q;ml=1-938 with initial composition X=0.7071 and Z=0. 01885. Further 
details of the physics employed by DSEP can be found in the remainder of this section and 
the two papers referenced in the first paragraph of this section paragraph. 
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Table 1. Solar Parameters 



Adopted Value Solar Model Value 



Age = 4.57 Gyr 
M = 1.9891 xlO 33 g 

R = 6.9598 xlO 10 cm Log(R/R )=O.OOOl 

L = 3.8418 xlO 33 erg/s Log(L/L )=O.OOO5 

R CZ /R Q = 0.713 ± 0.001 0.716 

Z/X = 0.0229 0.02288 



References. — Bahcall et al. (2005); Grevesse & 
Sauval (1998) 
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Table 3. Minimum Mass in M Q for Convective Core Overshoot 



Y=0.245+1.54Z Y=0.33 Y=0.4 



[Fc/H] [a/Fe}= -0.2 0.2 0.4 0.6 0.8 0.4 0.4 



-2.5 1.6 1.6 1.6 1.5 1.3 1.3 1.5 1.5 1.5 1.5 

-2.0 1.4 1.4 1.4 1.4 1.3 1.3 1.5 1.5 1.5 1.5 

-1.5 1.4 1.4 1.4 1.4 1.3 1.3 1.4 1.4 1.4 1.4 

-1.0 1.4 1.4 1.4 1.4 1.2 1.2 1.4 1.4 1.4 1.4 

-0.5 1.2 1.2 1.2 1.2 1.2 1.0 1.2 1.2 1.2 1.2 

0.0 1.2 1.2 1.1 1.1 1.0 0.9 1.2 1.1 1.2 1.1 
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The adopted convective core overshoot (CCO) treatment is that of Demarque et al. 
(2004). The basic idea behind the method is that the extent of CCO depends on the size 
of the convective core. The lowest mass main sequence stars that develop convective cores 
have relatively small ones. As stellar mass increases beyond this minimum, the convective 
cores become fully established and the amount of overshoot increases. 

The extent of CCO is treated as a fraction of the pressure scale height at the core 
boundary. At the minimum mass (listed for all compositions in Table 3) the amount of CCO 
is 5% of the pressure scale height. Between the minimum and 0.2 M above the minimum, 
the extent of CCO is ramped linearly from 5% to 20%. For masses greater than 0.2 M Q 
above the minimum, the extent of CCO is held constant at 20% of the pressure scale height. 
CCO plays a role in only the youngest isochrones presented in this paper, those with ages 
less than about 5 Gyr. 

2.2. Microscopic Diffusion and Gravitational Settling 

Chaboyer et al. (2001) described the implementation of diffusion and gravitational set- 
tling in DSEP; the following is a brief summary. The basic formalism is that of Thoul et al. 
(1994) but diffusion has been inhibited in the outer 0.1 M Q . In the outermost 0.05 M diffu- 
sion is completely stopped and in the 0.05 M below that diffusion is linearly ramped from 
zero to the full effect. This approach is meant to satisfy the lack of observational evidence 
for diffusion in metal poor stars (Gratton et al. 2001; Ramirez & Cohen 2003, e.g.), though 
recently evidence for some diffusion has been reported in NGC 6397 (Korn et al. 2006). 

2.3. Equation of State 

The basic equation of state (EOS) employed by DSEP is a general ideal gas equation 
of state with the Debye-Huckel correction and performs sufficiently well for stellar models 
with masses greater than about 0.7 M (Chaboyer & Kim 1995). Because this EOS neglects 
many non-ideal effects that become important at lower temperatures and higher densities 
(such as are found in low mass main sequence stars) the basic EOS was used only for models 
with masses greater than 0.8 M . 

For stellar models with 0.1 < M < 0.8 M the detailed EOS code developed by Irwin 
(2004), FreeEOS (in the EOS4 configuration), was used. FreeEOS takes into account many 
non-ideal effects and is generally superior for high density, low temperature environments. 
Furthermore, FreeEOS takes into account the makeup of Z. Stellar models produced with this 
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EOS are in good agreement with observations of low mass stars (see §5). Irwin has written 
extensively about the numerical methods and physical models employed by FreeEOS; the 
interested reader should consult the FreeEOS website 2 . 

Beginning at approximately 0.7 M , stellar evolution tracks produced with FreeEOS 
and the default EOS begin to differ in the H-R diagram. As stellar mass decreases, these 
differences increase, becoming substantial at 0.1 M . The current version FreeEOS (2.0) 
does not operate in the parameter range occupied by main sequence models below 0.1 M Q . 

2.4. Opacities 

Opacities were adopted from two sources 3 . In both cases, separate sets of tables were 
generated for each of the [cc/Fe] mixtures listed above. The OPAL opacities (Iglesias & 
Rogers 1996) were used at high temperature (log T > 4.5) while opacities computed for this 
effort following Ferguson et al. (2005) were used at low temperature (log T < 4.3). Between 
log T = 4.3 and 4.5 the opacity was ramped between the two sources. Ferguson et al. (2005) 
have shown that OPAL opacities agree with their own to within 5% in this temperature 
range. 

During core and shell He fusion (HB and AGB evolution) the OPAL Type 2 opacity 
tables were used. The Type 2 opacity tables allow mass fractions of two individual elements 
(in this case C and O) to vary independent of overall Z. These opacity tables are particularly 
helpful in the case of He-burning evolution where the He core gradually becomes a mixture 
of almost entirely C and O over the course of ~100 Myr. 

Conductive opacities were adopted from Hubbard & Lampe (1969) in the non-relativistic 
limit and Canuto (1970) in the relativistic limit via the fitting formulas of Sweigart (1973). 

2.5. Surface Boundary Conditions 

It is typical in stellar evolution codes to see surface boundary conditions treated as 
simple mathematical models of surface temperature as a function of optical depth. The 
almost universal choices are the Eddington approximation or the empirical fit to the Sun 



2 http://freeeos. sourceforge.net 

3 OPAL opacities — http:/ /www-phys. llnl.gov/Research/OPAL/ 
Ferguson et al. (2005) opacities — http://webs.wichita.edu/physics/opacity/ 
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done by Krishna Swamy (1966). 

However, particularly in the case of cool, low mass main sequence stars it is common 
practice to derive surface boundary conditions from model atmosphere calculations (Baraffe 
et al. 1997). Model atmosphere boundary conditions are an important ingredient in models 
of cool or very low mass stars (M < 0.4 M Q ) but make little or no difference for hotter, more 
massive stars. 

DSEP has been modified to accept surface boundary conditions derived from PHOENIX 
model atmospheres (Hauschildt et al. 1999a,b). One of the strengths of the PHOENIX code 
is attention to the effects of molecules and grains that appear at low temperature. Indeed, the 
low temperature opacities used in this paper were computed with a modified version of the 
PHOENIX code (Ferguson et al. 2005). A grid of stellar atmosphere models was produced 
covering the range of [Fe/H] and [a/Fe] listed above based on the Grevesse & Sauval (1998) 
solar abundances. The grid spans surface gravity -0.5 < log g < 5.5 and effective temperature 
2,500 K < T cff < 10,000 K. For the hottest tracks, model atmosphere grids from Castelli & 
Kurucz (2003) were supplemented for T eff > 10,000 K. A coarser grid of model atmospheres 
was generated to explore the effect of enhancing the He abundance on atmosphere structure 
and synthetic flux. While additional He in the atmosphere has a measurable effect on the 
gas pressure there is essentially no change in the synthetic flux. Because He is a noble gas, 
it is difficult to excite for T e g below 10,000 K and thus contributes little to the flux. At 
present, He-enhancement is treated as an adjustment to the gas pressure that serves as the 
surface boundary condition but is ignored in the synthetic colors. 

The surface boundary condition is satisfied by specifying pressure as a function of T c g- 
and log g. To obtain a pressure from each model atmosphere, the location within the 
atmosphere structure where T=T c g was found and the gas pressure extracted from that 
location. This approach differs from those in which the pressure is extracted from the model 
atmosphere at constant optical depth. Before deciding to follow the adopted method, tests 
were performed by comparing surface boundary condtions extracted from the atmospheres 
at constant optical depth (r=2/3 and 10) but the tests showed that the T=T e fr method 
produced the best results when compared with observations. The physical approximations 
made by model atmosphere and stellar evolution codes are different enough than any choice 
of where to connect the interior to the atmosphere will encounter some inconsistencies. 

Calibrated solar models using the Eddington, Krishna Swamy, and PHOENIX surface 
boundary conditions were produced in order to compare the three boundary conditions in 
the H-R diagram. The Krishna Swamy model required the largest mixing length (q;ml=2.02) 
and has the hottest red giant branch. The Eddington model required the smallest mixing 
length (o:ml=1-73) and has the coolest red giant branch. The PHOENIX model was located 
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Effective Temperature (K) 



Fig. 1. — The evolution of calibrated solar models from the pre- main sequence to the tip 
of the red giant branch the in the theoretical H-R diagram (except for K-S where no pre- 
main sequence is plotted). Each model employs a different surface boundary condition 
and requires a different solar-calibrated mixing length. The PHOENIX model atmosphere 
boundary condition adopted in this paper lies between the Eddington and Krishna Swamy 
models both in terms of the mixing length and, as a result, the temperature of the red giant 
branch. Note also that while the Eddington and Krishna Swamy models have a similar 
amount of curvature above log L/L =l the PHOENIX model red giant branch is relatively 
straight. 
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in between the others (oml=1-94) with a slightly flatter red giant branch in the H-R diagram 
(see Figure 1). 

2.6. Nuclear Reactions and Neutrino Cooling Rates 

Nuclear reaction rates were adopted from (Adelberger et al. 2002) with two exceptions: 
first, the 14 N(p,7) 15 reaction (Imbriani et al. 2004) and second, the 12 C(a,7) 16 reaction 
(Kunz et al. 2002). Neutrino cooling rates were adopted from Haft et al. (1994). 

2.7. He-Burning Evolution 

Stellar evolution calculations that follow the He flash in low mass stars are possible 
but numerically difficult and time consuming. Furthermore, studies have shown that so- 
called pseudo-ZAHB (zero age horizontal branch) models that are concocted from red giant 
branch tip models are nearly indistinguishable from full He flash calculations in terms of the 
subsequent evolution (Vandenberg et al. 2000; Piersanti et al. 2004; Serenelli & Weiss 2005). 
Beginning with a number of these pseudo-ZAHB models a complete set of core and shell 
He-burning tracks were created based on the red giant branch tip surface compositions listed 
in Table 4 and He core masses listed in Table 5. Since He core mass varies slightly with age 
or stellar mass, the core mass was taken from the stellar track that reached the red giant 
branch tip closest to 12 Gyr in each case, if no track for a given composition reached the 
tip within 11-13 Gyr the surface composition and He core mass were interpolated between 
the two nearest tracks. The core and shell He-burning tracks allow for the population of 
the horizontal branch and asymptotic giant branch through the use of individual tracks or 
the construction of synthetic horizontal branch models (for more on this see §4.4). The 
remainder of this section describes the physics in DSEP specific to core and shell He-burning 
processes. 

Semiconvection has long been known to play an important role in stars with steep 
composition gradients such as occur in stars undergoing core He fusion. Instead of including 
a full treatment of a semiconvection zone, DSEP uses a simple model to treat the expansion 
of the core boundary during He-burning developed by Castellani et al. (1971). Thermal 
energies have been artificially reduced when the core He mass fraction falls below 0.1 to 
discourage small loops in the H-R diagram known as breathing pulses as suggested by A. 
Sweigart (private communication) and first noted by Dorman & Rood (1993). Breathing 
pulses occur when the convective core expands outward into He-rich material and are widely 
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^Y init and Z init refer to the initial Y and Z values, for the initial compositions see Table 2. The Z values in this table 
differ from those in Table 2 as a result of diffusion. 
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believed to be unphysical artifacts of the models. Neglecting thermal energies in the core 
during the final 10% of core He fusion effectively freezes the size of the core thereby removing 
breathing pulses. 

High temperature opacities that allow for varying levels of C and O (as described in 
§2.4) allow for correct opacity coverage of stellar cores during the entire core and shell He 
fusion phases. 

Time steps during core He fusion are based on the depletion of He in the core such 
that the mass fraction Y core cannot fall by more than a fixed amount during each time step. 
When Ycore falls to 0.1 then the time step is set by the percentage change in Y core rather 
than by a a fixed amount. During shell He fusion (the AGB phase) time steps are further 
constrained by the total He-burning luminosity: the sum of the triple-a and 12 C(o;,7) 16 
luminosities. These time steps are used in addition to the time steps that govern the H-shell. 

3. Color-T cff Transformations 

The stellar evolution models described in the last section have been transformed to the 
observational plane by two different methods. One is based on the semi-empirical colors 
and bolometric corrections of Vandenberg & Clem (2003) and the transformation equations 
found in Appendix D of Sirianni et al. (2005), specifically those in Tables 22 and 24, hereafter 
referred to as the Empirical transformation. The other is based entirely on synthetic fluxes 
from PHOENIX (the same model atmosphere grids described in §2.5) with the appropriate 
filter transmission curves and is hereafter referred to as the Synthetic transformation. 

Both transformations have been normalized to Vega and include ground-based B, V and 
Table 5. He Core Mass in M for Zero Age Horizontal Branch Models 



Y=0.245+1.54Z Y=0.33 Y=0.4 



[Fe/H] [a/Fe]= -0.2 0.2 0.4 0.6 0.8 0.4 0.4 



-2.5 0.506 0.505 0.504 0.502 0.500 0.497 0.489 0.485 0.474 0.473 

-2.0 0.500 0.499 0.498 0.494 0.492 0.491 0.483 0.479 0.469 0.468 

-1.5 0.495 0.493 0.491 0.489 0.488 0.486 0.478 0.474 0.464 0.463 

-1.0 0.489 0.488 0.486 0.484 0.481 0.478 0.473 0.469 0.460 0.458 

-0.5 0.483 0.482 0.480 0.477 0.473 0.469 0.469 0.465 0.457 0.455 

0.0 0.475 0.473 0.469 0.464 0.458 0.447 0.465 0.460 0.454 0.451 
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I magnitudes on the Bessell (1990) system and the HST filters F606W and F814W on both 
the WFPC2 and ACS-WFC systems from Sirianni et al. (2005) 4 . 

3.1. Empirical Transformation 

Vandenberg & Clem (2003) introduced a new set of semi-empirical color transformations 
and bolometric corrections for standard B, V, R, and I bandpasses. The colors and bolometric 
corrections were empirically constrained so that the isochrones of Berbusch & Vandenberg 
(2001) were able to accurately reproduce the color-magnitude diagrams (CMD's) of various 
Galactic globular, open clusters, and assorted other stars spanning the range of metallicities 
commonly observed in the Milky Way. The BCy, B-V, and V-I of Vandenberg & Clem 
(2003) make up the ground-based portion of the Empirical transformation. 

As part of an exhaustive study of the ACS, Sirianni et al. (2005) derived equations to 
convert between ground-based BVRI magnitudes and HST ACS-WFC (and HRC) filters as 
well as between ACS and WFPC2 magnitudes. The V and I magnitudes from Vandenberg 
& Clem (2003) were used to obtain ACS F606W and F814W magnitudes using the (V,V- 
I)^F606W and (I,V-I)^F814W equations. The ACS magnitudes were then used to obtain 
the corresponding WFPC2 magnitudes. The VegaMAG zeropoints were used in each case. 

3.2. Synthetic Transformation 

Filter transmission curves in B, V, and I from Bessell (1990) were used along with curves 
from both ACS-WFC and WFPC2 filters F606W and F814W from Sirianni et al. (2005) to 
create a synthetic color-T eff transformation. The synthetic colors were generated with fluxes 
from the same PHOENIX model atmospheres used for the surface boundary conditions 
described in §2.5. These colors account for the effects of a-enhancement — something not 
possible with the empirical transformation. The effects are small but noticeable, particularly 
below ~4,000 K. 

As mentioned in §2.5, model atmospheres with enhanced He were produced to properly 
treat He-enhancement in the surface boundary condition. Tests have been performed to ad- 
dress the possible changes to synthetic colors due to the enhancement of He in the atmosphere 
for the Y=0.33 and 0.40 isochrones. As discussed in §2.5, the flux is virtually unchanged by 
enhancing He and this carries through to the colors. Therefore He-enhancement has been 



Readers interested in obtaining isochrones in other filters should contact the authors. 
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omitted from the Synthetic color transformation. 

The procedure followed to derive the Synthetic transformation closely follows that of 
Bedin et al. (2005). The bolometric correction in V and all colors were normalized by 
assuming Vega has V=0.026 and all colors equal to zero. 



Isochrones constructed from the two color transformations tend to agree within ~0.05 
mag in all filters (the only noteworthy exception occurs at high Z where the Synthetic B is 
flawed). Figure 2 compares the two color transformations over the full range of the isochrones 
for the case of [Fe/H]=-2.0, [a/Fe]=+0.4 at 12 Gyr in the left panel and [Fe/H]=0, [a/Fe]=0.2 
at 5 Gyr in the right panel. In each panel, the transformed isochrones were plotted with no 
adjustments in either direction, thus any difference is due entirely to the color transformation. 
The Empirical isochrones (dashed line) generally have a narrower color difference between 
the turn off and the base of the red giant branch than the Synthetic isochrones (solid line). 
The red giant branches move closer together until they overlap (or nearly so) at the tip. 
There is good agreement for the first ~5 mag below the turn off and then the Empirical 
transformation moves to the red of the Synthetic. 

Further side-by-side comparisons of the two color transformations for a variety of metal- 
licities and ages, including comparisons to photometry with distance and reddening esti- 
mates, can be found in Sarajedini et al. (2007). 



The following section lists and describes the different types of models accompanying 
this paper, either as data files or files that may be generated with one of several computer 
programs distributed along with the data files. The contents of each type of data file and 
the use of the accompanying programs are described in the Appendix. 



3.3. 



Discussion 



4. 



Stellar Evolution Models 



4.1. 



Stellar Evolution Tracks 



Stellar evolution tracks of individual masses for each composition are available in two 
formats: one listing the theoretical quantities and one listing basic global parameters and 
absolute magnitudes in each of the seven filters listed in the previous section. The tracks 
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F606W-F814W F606W-F814W 



Fig. 2. — Color comparison for the ACS bands: Left panel — Isochrones with [Fe/H]=-2.0, 
[ a /Fe]=+0.4 at 12 Gyr. Right panel— Isochrones with [Fe/H]=0, [a/Fe]=+0.2 at 5 Gyr. The 
general trend is for the Empirical colors to be bluer on the giant branch and redder on the 
lower main sequence than the Synthetic colors. (Empirical — dashed line; Synthetic — solid 
line) 
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Fig. 3. — He-burning tracks for [Fe/H] =-2 (top), -1 (middle), and (bottom) all with scaled- 
solar abundances. All masses from the minimum on the horizontal branch (see Table 5 for 
He core masses) of 0.51 M at [Fe/H]=-2, 0.5 M at [Fe/H]=-1, and 0.48 M at [Fe/H]=0 
up to 0.9 M are shown. The left side shows the Log T e g vs. Log L/L Q and the right side 
the V vs. V-I color- magnitude diagram using Empirical color transformation. The [Fe/H]=0 
color-magnitude diagram has been truncated in order to show the hotter stars more clearly, 
the tracks extend out to V-I~4. 
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include evolution from the pre-main sequence to the tip of the red giant branch from 0.1 M Q 
< M < 1.8 M Q in steps of 0.05 M . 

Separate tracks follow the He-burning evolution from the ZAHB to the onset of the 
thermal pulsation (TP-AGB) phase; these tracks are produced at masses from ~0.01 M Q 
above the He core mass (see Table 5) up to 0.51 M in steps of 0.01 M , from 0.52 to 0.56 
M in steps of 0.02 M , from 0.6 to 0.9 M in steps of 0.1 M , and finally from 1 to 1.75 M 
in steps of 0.25 M . Figure 3 shows all of the He-burning tracks from the minimum mass 
up to 0.9 M for [Fe/H]=-2, -1, and with scaled-solar abundances both in the theoretical 
H-R diagram and the V vs. V-I color- magnitude diagram. The [Fe/H]=0 tracks have been 
truncated in the color-magnitude diagram in Figure 3 in order to more clearly show the 
hotter stars, these tracks extend out to V-I~4. 

Each track file containing theoretical quantities lists for each time step: age in years, log 
effective temperature in Kelvin, log surface gravity (CGS), log luminosity in L , log radius 
in R , core He mass fraction Y core , core heavy element mass fraction Z core , surface heavy 
element to hydrogen ratio (Z/X) SMr j, H luminosity in L , He luminosity in L , He core mass 
in M , and CO core mass in M . All logarithms are base 10; see Table 1 for the adopted 
solar values. 

Each track file containing observable quantities lists at each time step the first four 
quantities from the theoretical track file (age, log T eff , log g, and log L/L ) plus absolute 
magnitudes in each of the seven filters listed above. Samples of both types of track file are 
given in the Appendix. 

4.2. Isochrones 

Isochrones were built using the equivalent evolutionary phase (EEP) principle. Each 
isochrone consists of ~230 points extending from the lowest mass point at 0.1 M up to the 
tip of the red giant branch. EEPs are assigned as a function of core He mass fraction (Y core ) 
on the main sequence and as a function of He core mass on the red giant branch. Isochrones 
are provided for ages ranging from 2 to 15 Gyr in increments of 0.5 Gyr for all compositions. 

Figure 4 demonstrates the effects of He-enhancement at [Fe/H]=-1, [a/Fe]=0. Isochrones 
are shown for 4, 8, and 12 Gyr with Y init ~0.25 (solid lines) and Y init =0A (dashed lines). 
Enhancing He at the expense of H decreases the opacity and increases the mean molecular 
weight. The end result is an isochrone with Y=0.40 is significantly hotter at a given age and 
Z than its counterpart with Y~0.25 as is demonstrated in Figure 4. 
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Effective Temperature (K) 

Fig. 4. — The effect of He-enhancement at constant [Fe/H]=-1 and [a/Fe]=0 in the H-R 
diagram. Ages shown are 4, 8, and 12 Gyr. The solid lines have Y ini4 =0.248 while the 
dashed lines have Y init =0.4. 



-22 - 



The model grid has a relatively fine spacing in both age and [a/Fe], such that parameters 
of clusters being analyzed can be well constrained by the ages and [a/Fe] values provided. 
The only parameter that could be more tightly constrained by observations is [Fe/H] and 
so a computer program has been distributed with the isochrones to allow interpolation in 
[Fe/H] (see the Appendix for details). The interpolation code has been extensively tested 
and provides results that are reasonably accurate over the range of metallicities provided in 
the isochrone database. 



4.3. Luminosity Functions 

Rather than adopt one or more standard filters and initial mass functions (IMF) and 
provide tabulated luminosity functions the distribution of isochrones comes with a program 
to generate luminosity functions from any of the isochrone files. Luminosity functions may 
be constructed from any bandpass with a choice of bin size and IMF exponent. The IMF 
is assumed to be a simple power law where dN/dM oc M x . Figure 5 demonstrates the 
luminosity functions when compared to the observed luminosity function of M92 (Paust et 
al. 2007). The isochrones on which the luminosity functions were based have [Fe/H] =-2.3 
and [a/Fe]=0.4 at ages 12, 13, and 14 Gyr. The assumed IMF is the Salpeter IMF (x=- 
2.35), a good fit in this case. The distance modulus is DMy=14.6 as used in Figure 11 for 
the Synthetic isochrone. The theoretical luminosity functions extend approximately seven 
magnitudes fainter than the data in Figure 11. 

4.4. Synthetic HB Models 

The formalism for constructing synthetic horizontal branch (SHB) models is based on 
that of Lee, Demarque, & Zinn (1990). Accompanying the He-burning evolutionary tracks is 
a program for creating SHB models. The SHB models are constructed from the tracks based 
on a set of parameters supplied by the user. These parameters include the composition 
([Fe/H], [a/Fe] ,Y), average HB mass (< M H b >), the mass of the red giant branch tip 
(Mrgbt, effectively the upper bound to the mass distribution), the mass dispersion on the 
HB (a M ), and the total number of stars (N). Based on these parameters, the program 
constructs a random Gaussian distribution of stellar masses, chooses a random age, and 
interpolates within the tracks for mass, age, and composition. Interpolation is linear within 
the tracks (age) and quadratic between the tracks (mass and composition). Because the 
actual mass distribution is finite and truncated on both ends by the finite mass range the 
code calculates the actual values of <M HB > and (The and reports these at the end. 
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v 

Fig. 5. — Example of theoretical luminosity functions compared to the observed luminosity 
function of M 92 (Paust et al. 2007). The theoretical luminosity functions shown were derived 
from Synthetic color isochrones at [Fe/H]=-2.3 and [a/Fe]=0.4 for ages 12, 13, and 14 Gyr 
with a Salpeter IMF (see also Figure 11). The assumed distance is DMy=14.6. The full 
theoretical luminosity functions extend several magnitudes fainter than the data shown in 
the plot. 
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Fig. 6.— SHB models with N=1000 stars for [Fe/H]=-2.3 (top) and -0.75 (bottom). The top 
panel represents a blue HB in a metal poor cluster while the bottom panel represents a red 
HB in a metal rich cluster. The left hand side shows the SHB models in the color-magnitude 
diagram while the right hand side shows the corresponding luminosity function. 
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In addition to the SHB generating program, another program is supplied to generate 
luminosity functions and fiducial lines from each individual SHB model. The program works 
on the raw output from the SHB code by binning the SHB model data, determining the 
mean and standard deviation of magnitude and color in a given bin, and writing results from 
all bins to a file. Use of this program in conjunction with the isochrone-based luminosity 
function code allows the user to create full luminosity functions that include all phases of 
the evolution covered by the models presented in this paper. 

Figure 6 demonstrates the effect of metallicity and the mass variables mentioned above 
in the color- magnitude diagram for two SHB models on the right. The left side of Figure 
6 shows luminosity functions constructed from the corresponding SHB model. The upper 
panels show the CMD and number distribution of a blue, metal poor SHB model with 
[Fe/H]=-2.3. The lower panels show the same information for a red, metal rich SHB model 
with [Fe/H]=-0.75. These models, shown alone here, are the same used in Figures 11 and 
12 below. 



5. Comparisons 

5.1. Comparisons with other isochrones 

Several groups currently produce stellar evolution tracks and isochrones similar to those 
presented in this paper. Here the new DSEP isochrones are compared to isochrones from the 
other groups in both the theoretical H-R diagram and in the V vs. V-I color-magnitude dia- 
gram. Comparisons in the H-R diagram allow for differences in the luminosity and effective 
temperature scales to be assessed while differences in the color-magnitude diagram display 
additional differences due to the color-T cff transformations adopted by each group. The 
DSEP isochrones shown in the following figures employ the Synthetic color transformation. 
Figure 2 shows the differences between Synthetic and Empirical colors. 

Comparisons were performed between DSEP and BaSTI, Padova, Victoria-Regina, and 
Y 2 at 10 Gyr with Z=0.0001 assuming a scaled-solar heavy element mixture. Scaled-solar 
was chosen because the definition and amount of a-enhancement is not consistent amongst 
the groups. Neither the primordial He abundance nor the adopted AY/AZ is consistent, at 
Z=0.0001: BaSTI has Y=0.245, DSEP has Y=0.02452, Padova has Y=0.23, Victoria-Regina 
has Y=0.2352, and Y 2 has Y=0.2302. Aside from differences in composition, the groups do 
not all use the same assumptions concerning the surface boundary condition, solar-calibrated 
mixing length, equation of state, gravitational settling, mass loss, and other physics: each of 
these discrepancies will result in differences amongst the isochrones at constant age and Z. 
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Figure 7 compares DSEP to BaSTI, Figure 8 to Padova, Figure 9 to Victoria-Regina, 
and Figure 10 to Y 2 . Both BaSTI and Padova include core He-burning evolution in their 
isochrones but these portions have been removed from the figures for clarity. 

In the H-R diagram, the DSEP isochrone is slightly hotter and more luminous than 
the others around the main sequence turn off and sub-giant branch but this does not carry 
over to the color-magnitude diagram in all cases due to the differences in adopted color- T ei j 
transformations. There is good agreement on the red giant branch tip location in the H-R 
diagram except with Padova which is hotter than all the others. 
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5.2. Comparisons with Data 

ACS photometry of two classic globular clusters, M92 and 47 Tuc, from Sarajedini et 
al. (2007) and Anderson et al. (2007, in preparation) that span a wide range of [Fe/H] have 
been chosen for comparison. Many of the brightest giants in Figures 11 and 12 include at 
least one saturated pixel. For these stars, the point spread function photometry used for 
unsaturated stars has been replaced by aperture photometry with excellent results. Full 
details of the photometry will be given by Anderson et al. (2007, in preparation). 

Isochrones for reasonable assumptions of the age and metallicity of M 92 and 47 Tuc 
are shown in Figures 11 and 12, respectively. The purpose of these plots is not to constrain 
any of these parameters for either cluster but merely to demonstrate the performance of the 
isochrones when compared with high quality photometry. To show the complete evolution- 
ary picture, ridge lines from the SHB models shown in Figure 6 are plotted alongside the 
photometry: the metal poor SHB model ([Fe/H]=-2.3) with M 92 and the metal rich SHB 
model ([Fe/H] =-0.75) with 47 Tuc assuming the adopted reddening and distance modulus 
in each figure. 

Figure 11 shows little difference between the color transformations and both provide 
reasonable fits over the extent of the data. The largest difference occurs on the subgiant 
branch: in this region the Empirical color transformation is clearly a better fit to the data. 
The behavior of the Synthetic transformation is common to most color transformations 
derived synthetic spectra. 

Figure 12 shows more clearly the difference between the two color transformations. The 
Empirical transformation provides a good fit from the red giant branch to approximately 3 
mag below the main sequence turn off. Below that point, the Empirical isochrone becomes 
increasingly redder than the data. In contrast, the Synthetic isochrone becomes bluer than 
much of the data a few magnitudes below the turn off. The difference between the Synthetic 
isochrone and the data is smaller than in the case of the Empirical isochrone, and in the 
opposite sense, but still visible. For this reason, the Synthetic color transformation is prefer- 
able for examining the lower main sequence, particularly at higher metallicity. Disagreement 
between the lower main sequence data and the isochrones is likely due to errors in both the 
effective temperature scale of the stellar evolution models and in the color transformations 
but it is not possible to quantify these errors. 




Fig. 7. — BaSTI (dashed line) and DSEP (solid line, with Synthetic color transformation) 
isochrones at 10 Gyr with Z=0.0001 (scaled-solar mixture). 
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Fig. 9. — Same as Figure 7 but comparing Victoria-Regina and DSEP isochrones. 
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1 2 1 2 

F606W-F814W F606W-F814W 



Fig. 11. — ACS data from M92 compared to isochrones with both the Empirical (left panel) 
and Synthetic (right panel) color transformations. Details are listed on each panel. Data 
from Sarajedini et al. (2007). The fiducial line from the metal poor synthetic horizontal 
branch model of §4.4, Figure 6, is plotted alongside both isochrones. 
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Fig. 12. — ACS data from 47 Tuc compared to isochrones with both the Empirical (left 
panel) and Synthetic (right panel) color transformations. Details are listed on each panel. 
Data from Sarajedini et al. (2007). The fiducial line from the metal rich synthetic horizontal 
branch model of §4.4, Figure 6, is plotted alongside both isochrones. 
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6. Conclusions 

The stellar evolution code DSEP has been modified to significantly improve the accuracy 
of isochrones on the lower main sequence and to self-consistently treat the outer layers of 
stars for a wide range of compositions. The synthetic and empirical color transformations 
employed have strengths and weaknesses and highlight the amount of uncertainty in current 
color-T c flf transformations in the color-magnitude diagram. 

A new set of stellar evolution models have been computed with DSEP for a range of 
initial He abundance, [Fe/H], and [ce/Fe] encompassing ages and abundances observed in 
Galactic globular clusters. The models are presented as isochrones (2-15 Gyr) and stellar 
evolution tracks (0.1-1.8 M ) with accompanying programs to create luminosity functions 
and SHB models suitable for analysis of old stellar populations. Details of the stellar evolu- 
tion calculations and all steps leading up to the final products have been discussed. These 
models should prove useful in analyzing the ACS Galactic Globular Cluster Survey data and 
in many other applications. The models extend to 0.1 M on the main sequence, making 
them useful for deep photometric studies such as the ACS Galactic Globular Cluster Survey 
and at the same time provide nearly full evolutionary coverage of the stars found in globular 
clusters and old open clusters. 
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A. Composition, Color Transformations, and File Names 

Composition is specified in the file names in terms of [Fe/H] as f eh and [a/Fe] as af e. 
The values of [Fe/H] and [a/Fe] are listed using the letters m for minus (<0) and p for plus 
(>0) followed by two digits for [Fe/H] and one digit for [a/Fe]. For example, [Fe/H] =-2.5, 
[a/Fe]=+0.4 would be represented by the string fehm25afep4. Zero is set to positive and 
thus [Fe/H]=0, [a/Fe]=0 is listed as f ehpOOaf epO. Tracks and isochrones with enhanced He, 
either Y=0.33 or Y=0.4, have an additional string following f ehXXXaf eXX in their filenames: 
y33 or y40, respectively. 

Color transformations are listed as a suffix to isochrone and evolutionary track file 
names. The synthetic color transformation is indicated by the suffix .phx and the semi- 
empirical color transformation by . cmd. As a final example, the synthetic color isochone 
with [Fe/H]=-1.5, [a/Fe]=+0.4, and Y=0.33 is called f ehml5af ep4y33.phx. 

B. Stellar Evolution Track Files 

The stellar evolution tracks begin on the pre-main sequence and terminate at the He 
flash or begin on the ZAHB and terminate at the onset of thermal pulsations on the AGB 
(TP-AGB). Tracks are provided for masses ranging from 0.1 to 1.8 M in increments of 
0.05 M . For the mass range that experiences the He flash, separate He-burning tracks are 
provided. The content of these tracks are identical to the H-burning tracks as described 
below. 

The track files come in two formats: one containing theoretical quantities and one 
containing observable quantities. The theoretical track files contain information about the 
initial composition and physics used in a header followed by data lines. Each data line lists 
age in years, log T cff , log g, log L/L Q , log R/R©, Y core , Z core , (Z/X) surf , L H /L Q , L He /L Q , 
Mne-core, and Mco-core- The quantity L# is the total H-fusion luminosity combining all 
branches of the p-p chain and CNO cycle luminosities while L# e combines luminosities from 
the triple-a and 12 C(a,7) 16 reactions. 

The naming convention for the stellar evolution tracks is to use the prefix m (for mass) 
followed by three integers representing the stellar mass (1 M Q — >■ mlOO) followed by com- 
position information and either the suffix . trk for theoretical quantities or the appropriate 
suffix for a color transformation specified above. The He-burning tracks use the same naming 
convention except the m has been replaced with hb. 

The header and first eight data lines of the 1 M Q , [Fe/H]=0, [a/Fe]=0, Y=0.27 stellar 
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evolution track are given below for both the theoretical file mlOOf ehpOOaf epO.trk in Table 
6 and the observable file with the PHOENIX color transformation mlOOf ehpOOaf epO .phx in 
Table 7. The columns labeled F606W and F8I4W refer to the ACS bands and the columns 
labeled F6062 and F8142 refer to the WFPC2 bands. 



C. Isochrones and Luminosity Functions 

Isochrone files begin with a header that lists composition and color transformation 
details as well as the number of ages in the file. The isochrone files distributed with this 
paper contain 27 ages from 2 to 15 Gyr in increments of 0.5 Gyr. The individual isochrones 
have their own headers consisting of two lines. The first line gives the age of the isochrone 
and the number of data points (EEPs) it contains. The second line lists the quantities 
provided on each line: EEP number, stellar mass, log T e ff, log g, log L/L Q , and the absolute 
magnitudes in all seven filters. The header and first six data lines of the [Fe/H]=0, [a/Fe]=0 
isochrone file fehp00afep0.cmd are printed in Table 8. The columns labeled F606W and 
F814W refer to the ACS bands and the columns labeled F6062 and F8142 refer to the 
WFPC2 bands. 

Output from the luminosity function program (isolf . f ) contains the same header as the 
isochrone file it was created from. Likewise, each luminosity function has a two line header 
that lists the age and number of data points followed by a description of the information 
printed in each line. The data lines list a bin number, the average absolute magnitude of the 
bin, the average color of the bin, the cumulative (logged) number of stars, and the (logged) 
number of stars in that bin. The cumulative star count begins with the brightest bin. The 
header and first six data lines of a luminosity function generated from the isochrone file 
(Table 8) are printed in Table 9. The luminosity function is created using the F606W (ACS) 
filter with color F606W-F814W, bin size of 0.1 mag, and Salpeter IMF. 

D. Synthetic Horizontal Branch Models and Luminosity Functions 

Output from the SHB code (shb.f) consists of a file containing one line for each star 
plus a short header. The header lists number of stars, [Fe/H], [a/Fe], and Yinit- The data 
lines include stellar mass, age, log T e s, log g, log L/L , and absolute magnitudes in all seven 
filters. Table 10 shows the header and first eight data lines from an SHB model with 1000 
stars at [Fe/H]=0, [a/Fe]=0, and Y=0.274. 

Output from shb . f may be processed by shblf . f to create a fiducial line and luminosity 



Table 6. Example of a Stellar Evolution Track with Theoretical Quantities 



M=1.00 X=7.0710E-01 Z= 


:1.8850E-02 COMP= 


: GS98 [a/Fe]= 


0.00 MIXL= 


:1.9380 CCO^O.OO 








Age (yrs) Log T 


Log g Log L 


Log R 


Y_core 


Z.core (Z/X)_surf 


L_H 


L_Hc 


M_He_core M_CO_core 



5.27000000E+02 
5.94200000E+02 
6.64760000E+02 
7.38848000E+02 
8.16640400E+02 
8.98322420E+02 
9.84088541E+02 
1.07414297E+03 



3.62873 
3.62878 
3.62884 
3.62890 
3.62897 
3.62904 
3.62911 
3.62919 



2.41053 
2.41191 
2.41335 
2.41487 
2.41645 
2.41810 
2.41983 
2.42163 



1.49574 
1.49456 
1.49335 
1.49209 
1.49077 
1.48940 
1.48796 
1.48646 



1.01367 
1.01298 
1.01225 
1.01150 
1.01071 
1.00988 
1.00902 
1.00811 



2.7402E-01 
2.7402E-01 
2.7402E-01 
2.7402E-01 
2.7402E-01 
2.7402E-01 
2.7402E-01 
2.7402E-01 



1.8850E-02 
1.8850E-02 
1.8850E-02 
1.8850E-02 
1.8850E-02 
1.8850E-02 
1.8850E-02 
1.8850E-02 



2.6658E-02 
2.6658E-02 
2.6658E-02 
2.6658E-02 
2.6658E-02 
2.6658E-02 
2.6658E-02 
2.6658E-02 



0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 



0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 



0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 



0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 
0.0000E+00 



I 

CO 

I 



Table 7. Example of a Stellar Evolution Track with Observable Quantities 



M/Mo= 


4.00 [Fe/H]= 0.07 


[a/Fe]= 


0.00 Y=0.274 








Age 


(yrs) Log T 


Log 


g Log L B 


V 


I 


F606W F814W F6062 F8142 



5.27000000E+02 
5.94200000E+02 
6.64760000E+02 
7.38848000E+02 
8.16640400E+02 
8.98322420E+02 
9.84088541E+02 
1.07414297E+03 



3.62873 2.41053 

3.62878 2.41191 

3.62884 2.41335 

3.62890 2.41487 

3.62897 2.41645 

3.62904 2.41810 

3.62911 2.41983 

3.62919 2.42163 



1.49574 3.1495 

1.49456 3.1517 

1.49335 3.1539 

1.49209 3.1562 

1.49077 3.1585 

1.48940 3.1609 

1.48796 3.1635 

1.48646 3.1662 



1.7757 0.3959 

1.7781 0.3988 

1.7806 0.4017 

1.7832 0.4047 

1.7859 0.4079 

1.7886 0.4112 

1.7916 0.4146 

1.7946 0.4182 



1.4607 


0.4112 


1.4632 


0.4140 


1.4657 


0.4169 


1.4684 


0.4199 


1.4712 


0.4231 


1.4740 


0.4263 


1.4771 


0.4298 


1.4802 


0.4334 



1.3968 0.4286 

1.3994 0.4315 

1.4019 0.4344 

1.4046 0.4374 

1.4074 0.4405 

1.4102 0.4438 

1.4133 0.4472 

1.4164 0.4508 



I 

CO 

oo 
I 



Table 8. Example of an Isochrone File 



NUMBER OF AGES=29 



MIX-LEN Y Z Zcff [Fe/H] [a/Fe] 

1.9380 0.2741 1.8850E-02 1.8850E-02 0.07 0.00 



PHOTOMETRIC SYSTEM: VfcC(2003) a +Sirianni et al.(2005) 



AGE= 2.000 EEPS=238 



EEP 


M/Mo 


LogTeff 


LogG 


LogL/Lo 


B 


V 


I 


F606W 


F814W 


F6062 


F8142 


4 


0.099162 


3.5048 


5.3118 


-2.9048 


16.5676 


14.8000 


11.7423 


14.1827 


11.6499 


14.0308 


11.6953 


5 


0.115037 


3.5048 


5.2473 


-2.7758 


16.2406 


14.4760 


11.4311 


13.8600 


11.3397 


13.7098 


11.3849 


6 


0.134056 


3.5048 


5.1810 


-2.6432 


15.9058 


14.1439 


11.1106 


13.5292 


11.0200 


13.3805 


11.0652 


7 


0.169755 


3.5076 


5.1016 


-2.4501 


15.3348 


13.5914 


10.6076 


12.9822 


10.5209 


12.8399 


10.5654 


8 


0.223896 


3.5162 


5.0395 


-2.2333 


14.5387 


12.8438 


9.9844 


12.2493 


9.9068 


12.1222 


9.9497 


9 


0.276573 


3.5236 


4.9959 


-2.0684 


13.9167 


12.2621 


9.5074 


11.6808 


9.4371 


11.5657 


9.4787 



a V&C2003 refers to Vandenberg & Clem (2003). 



Table 9. Example of a Luminosity Function File 



NUMBER OF AGES=29 



MIX-LEN Y Z Zeff [Fe/H] [a/Fe] 

1.9380 0.2741 1.8850E-02 1.8850E-02 0.07 0.00 



PHOTOMETRIC SYSTEM: VfcC(2003)+Sirianni et al.(2005) 



AGE= 
N 


2.000 BINS=158 IMF x= -2.35 
MF606W F606W-F814W 


LoglO(N) 


LoglO(dN) 


1 


-1.517 


1.7545 


0.1086 


0.1086 


2 


-1.417 


1.7064 


0.3138 


-0.1105 


3 


-1.317 


1.4011 


0.4090 


-0.2966 


4 


-1.217 


1.3380 


0.4921 


-0.2671 


5 


-1.117 


1.2829 


0.5655 


-0.2426 


6 


-1.017 


1.2381 


0.6276 


-0.2478 



Table 10. Example of a Synthetic Horizontal Branch File 



N= 1000 [Fe/H]= 0.00 


[a/Fe]=0.00 Y_i=0.274 
















M/Mo 


Age (yrs) 


Log Tcff 


Log g 


Log L/Lo 


B 


V 


I 


F606W 


F814W 


F6062 


F8142 


0.7771 


1.0639E+08 


3.6677 


2.2732 


1.6791 


2.1224 


0.9687 


-0.1202 


0.7106 


-0.1177 


0.7032 


-0.1000 


0.7920 


9.7970E+07 


3.6727 


2.3536 


1.6274 


2.1900 


1.0672 


0.0076 


0.8166 


0.0106 


0.8100 


0.0279 


0.8825 


6.9989E+07 


3.6760 


2.4404 


1.6004 


2.2174 


1.1160 


0.0741 


0.8700 


0.0774 


0.8639 


0.0944 


0.7812 


5.2351E+07 


3.6792 


2.4334 


1.5671 


2.2678 


1.1813 


0.1560 


0.9397 


0.1596 


0.9339 


0.1763 


0.8525 


1.0007E+08 


3.6694 


2.3350 


1.6638 


2.1377 


0.9968 


-0.0823 


0.7412 


-0.0796 


0.7341 


-0.0620 


0.8959 


1.1958E+08 


3.6137 


1.5453 


2.2529 


1.4486 


-0.0369 


-1.5517 


-0.3977 


-1.5599 


-0.4197 


-1.5359 


0.7885 


9.2760E+07 


3.6747 


2.3783 


1.6087 


2.2134 


1.1021 


0.0536 


0.8544 


0.0568 


0.8480 


0.0739 


0.8229 


1.0110E+08 


3.6695 


2.3239 


1.6602 


2.1459 


1.0050 


-0.0734 


0.7495 


-0.0707 


0.7424 


-0.0532 



-42- 



function from the SHB model. The shblf program creates luminosity functions and fiducials 
in any of the seven bands and with either a set number of bins or set bin size. It scans through 
and SHB model file and calculates the number of stars, the average and standard deviation 
of the absolute magnitude, and the average and standard deviation of the color in each bin 
and finally prints these data along with a short header. Table 11 show the header and first 
eight data lines of an SHB luminosity function derived from the SHB model described above 
(Table 10). 

E. Computer Programs and Summary 

Five computer programs, written in Fortran 77 are available along with the data files. 
These programs allow for the following manipulations of the data files: 

iso_interp_f eh.f : Isochrone interpolation in [Fe/H] at fixed [a/Fe], Y, and age 
iso_split.f: Create individual age isochrone files from one of the large, many-aged files 
isolf . f: Luminosity function construction from an isochrone file 
shb.f: Synthetic horizontal branch model construction with [Fe/H] interpolation 
shblf .f: Luminosity function/fiducial line construction from a SHB model file 

An accompanying README text file contains the text of this appendix. 

F. How to obtain the data files and programs 

The entire distribution of stellar evolution tracks, isochrones, and accompanying com- 
puter programs may be downloaded from http://stellar.dartmouth.edu/~models/ and the 
Multimission Archive at the Space Telescope Science Institute (MAST; http:/ /archive. stsci.edu/) 
as part of the ACS Galactic Globular Cluster Survey (GO-10775). 



Table 11. Example of a Synthetic Horizontal Branch Luminosity Function File 



N= 1000 [Fe/H]= 0.00 [a/Fe]=0.00 Y_i=0.274 
Bin Size= 0.0368 Number of Bins= 50 



N/Bin 


Bin Center 


<F606-F814> 


sigma 


<F606W> 


sigma 


5 


-0.82520300 


1.62138000 


0.11453257 


-0.82572000 


0.01088511 


2 


-0.78840900 


1.65370000 


0.11840000 


-0.79715000 


0.00715000 


4 


-0.75161500 


1.69665000 


0.24597608 


-0.74417500 


0.00733941 


3 


-0.71482100 


1.41780000 


0.02011765 


-0.71806667 


0.00549201 


2 


-0.67802700 


1.36825000 


0.03005000 


-0.68070000 


0.01010000 


2 


-0.64123300 


1.66125000 


0.34805000 


-0.64785000 


0.01055000 


3 


-0.60443900 


1.55713333 


0.35709536 


-0.60706667 


0.00978445 
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